Excitability in a nonlinear magnetoacoustic resonator 
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We report a nonlinear acoustic system displaying excitability. The considered system is a magne- 
tostrictive material where acoustic waves are parametrically generated. For a set of parameters, the 
system presents homoclinic and heteroclinic dynamics, whose boundaries define a excitability do- 
main. The excitable behaviour is characterized by analyzing the response of the system to different 
external stimuli. Single spiking and bursting regimes have been identified. 
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There is an increasing interest in the study of smart 
materials 0]. The term smart refers to a class of ma- 
terials that are highly responsive and have the inher- 
ent capability to sense and react according to changes 
in the environment. The most widely used classes of 
smart materials are ferroic, and include piezoelectrics, 
electrostrictors, magnetostrictors, and shape-memory al- 
loys. In this letter we consider a smart system based 
on the magnetostrictive effect. Magnetostriction is the 
material property that causes a material to change its 
dimensions under the action of a magnetic field or, in- 
versely, to generate a magnetic field when deformed by 
an external force. Magnetostrictive materials can thus 
be used for both sensing and actuation. 

One of the main goals of smart materials research is to 
design systems able to mimic biological behaviour. For 
example, artificial muscles have been recently developed 
based on electroactive polymers. Possibly the smartest 
biological system is the neuron, responsible among oth- 
ers of the signal processing in the brain. This process- 
ing capability relies to a large extend on a property on 
the individual neurons called excitability 0- The most 
characteristic feature of excitability is a highly nonlinear 
response (of all-or-nothing type) to an external stimu- 
lus: when the amplitude of the stimulus is below a given 
threshold, there is a weak reaction, while the response is 
strong, but independent of the strength of the stimula- 
tion, above the threshold. The interest in excitable sys- 
tems came first from chemistry (reaction-diffusion sys- 
tems) and biology (cardiac tissue and neural modeling) 
^ . More recently, excitability has been identified in cer- 
tain optical systems, as the CO2 laser with saturable ab- 
sorber 4] or the semiconductor laser with optical feed- 
back 5]. In this letter we report on the first, to the best 
of our knowledge, acoustic system displaying excitabil- 
ity. The system also displays a rich complex behaviour 
including an scenario of homoclinic/heteroclinic dynam- 
ics, which is on the basis of the excitable character of the 
system. 

The considered physical system consists in a magne- 
tostrictive material (e.g., a ferrite) in the form of a rod 
with transverse section S and length L, with parallel 
and flat lateral boundaries, which is driven by an ex- 
ternal magnetic field H^xt = Hq + Hq{t) parallel to the 



axis of the rod. The oscillating magnetic field Hq{t) = 
h co8{2ujt) induce a temporal modulation of the sound 
velocity in the material at frequency 2uj, which is respon- 
sible for the parametric excitation of an acoustic wave in 
the material at half of the driving frequency. A stand- 
ing wave is formed since the polished boundaries define 
a high-Q resonator for the sonic beam. This effect has 
been considered in a number of works 0, I3j where 
a radiofrequency magnetic field was used to excite para- 
metrically an ultrasonic wave. In these works, a special 
attention was paid to the phenomenon of phase conjuga- 
tion (PC) of an incident acoustic beam, with important 
applications in e.g. acoustic microscopy 0. 

In the simplest scheme, the oscillating magnetic field 
is produced by a coil with n turns surrounding the ma- 
terial. The coil provides the inductance of an electric 
RLC series circuit, driven by an ac source at frequency 
2uj and variable amplitude £. Owing to magnetoelas- 
tic coupling, elastic deformations in the magnetostric- 
tive material result in an additional magnetic field. Fol- 
lowing Isl 01) we consider that the dominant contribu- 
tion of the magnetoelastic interaction is quadratic in the 
particle displacements Hint{t) — —a (u{Y,tY^, where 
a — {2k)^{d\nv/dH) is the coupling coefficient propor- 
tional to the modulation depth of sound velocity |S| , and 
the brackets indicate a spatial average over the material 
volume. Taking into account the magnetoelastic contri- 
bution, the effective magnetic excitation in the material 
takes the form H = Hext{t)+Hint{t). The resulting mag- 
netic induction B = fi {H) H is in general a nonlinear 
relation, which to the leading order can be written as 
B = fiH + ^^qX^^^H^ 0, where /i is the linear perme- 
ability of the material and x'-^'' the third order magnetic 
susceptibility, which in turn depends on the frequency. 

A nonlinear equation for the electric circuit can be ob- 
tained under the previous assumptions, and neglecting 
nonresonant terms and those higher than quadratic in 
It reads 
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where q is the charge in the capacitor, related with 
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the current as / = dq/dt, £ is the coil inductance and 
Hg = nl. The last two terms result from the nonlinear- 
ities related to magnetoelastic interaction and magnetic 
nonlinearity. 

The acoustic field obeys the wave equation with a 
source (coupling) term proportional to the exciting mag- 
netic field y • In terms of the charge it reads 
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We consider solutions of Eqs. Q and (0) in the form of 
quasi-harmonic waves, i.e., whose amplitudes are slowly 
varying in time. In this case we can write 



(3a) 



u (r, i) = - [[/ (t) exp [iuit) + c.c] g (r^) sin [kz) , (3b) 

where u (r, t) is an eigenmode of the acoustical res- 
onator, k = m-K / L define the cavity resonances, and 
2a; = is the pumping frequency at the circuit 

resonance. In the slowly-varying envelope approxima- 
tion, the amplitudes obey the evolution equations (see 
for details) 



'^=r-X + Y^+ii^\Y\^X, 
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where X and Y are normalized values of the charge in 
the capacitor and the amplitude of the ultrasonic field, re- 
spectively, which relate to the physical variables through 
the transformations X = (v^na/2jjj) Q and Y = 
J (van) (iiL/C^ij^qY/^U. The pumping term V ~ 
{v'^na/Auj'yuR)£ is proportional to the voltage of the ac 
source. The parameter 7 = IuHq represents the ra- 
tio between losses, being 7q — R/2C and 7^ the electric 
and acoustic decay rates respectively. The last parameter 
is introduced phenomenologically, and take into account 
the losses due mainly to radiation from the boundaries. 
A dimensionless time r = ^qt has been also defined. 
Finally, rj = AljHq^ux''^^ joaP' pi^ remains as the single 
parameter of nonlinearity, with /x^ = n/ ijlq- Due to the 
number of parameters involved in 77, it can be varied over 
a wide range of values. Note that Eqs. |@J also possess 
the Z2 (refiection) symmetry {X,Y) — > [X^—Y). 

The stationary solutions of Eqs. I@J and their stability 
have been analysed in ■ This analysis showed that the 
trivial solution X = V, Y — 0, where the acoustic sub- 
harmonic field is switched-off, experiences a subcritical 
pitchfork bifurcation at P = 1, the parametric genera- 
tion threshold, and gives rise to a finite amplitude solu- 

tion \X\ = 1, = (1 ± ^7^2(1 + ,y2)_^2)/(i + 

Subcriticality implies a domain of coexistence between 
the trivial (rest) state and the finite-amplitude solution 



when Vt < V < 1, where Vt = v/ + ^he turn- 
ing point. For some values of the parameters, the finite 
amplitude solution can experience a Hopf (selfpulsing) 
instability, whose analytical expression was also given in 
111]. The boundaries of these local instabilities are de- 
picted in Fig.l. 
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FIG. 1: Bifurcation diagram of Eqs. (4) for 7 = 0.1. Lo- 
cal bifurcations are denoted by TP (fold or turning point), 
HB (Hopf bifurcation), and the line V = \ corresponding to 
the subcritical pitchfork bifurcation (parametric generation 
threshold). Three curves labeled as (a), (b) and (c) denote 
the locus of global bifurcations (see text for details). 

In ^3 it was also shown that the system possess ho- 
moclinic (global) dynamics at a particular value of the 
nonlinearity parameter r\. In this letter we perform a 
detailed two-parameter numerical analysis of Eqs. I@J. 
For this aim, we reduce the dimensionality of the pa- 
rameter space by fixing the value of the relative losses 
7 = 0.1, which is in correspondence with typical exper- 
imental conditions 0- The results for other values of 7 
are qualitatively similar. 

The bifurcation diagram shown in Fig. 1 reveals that 
the system presents a complex scenario of global dynam- 
ics. Besides the Hopf bifurcation leading to self-pulsing, 
three different global bifurcations have been identified 
in this system. These bifurcations can be detected nu- 
merically by computing the period of the limit cycles, 
since close to a homoclinic/heteroclinic point the pe- 
riod diverges to infinity as T = — ■iln(7'/^ — V), where 
Vfi — V measures the distance to the homoclinic bifurca- 
tion (which is assumed small) and A is the eigenvalue in 
the unstable direction of the saddle point • In Fig- 1 , 
infinite period (Tea) bifurcations are the curves labeled 
(a)-(c). Curve (a) corresponds to a gluing (double homo- 
clinic) bifurcation This bifurcation is characteristic 
of systems with Z2 symmetry and is mediated by a sad- 
dle point, which in this case corresponds to the trivial 
state. The gluing bifurcation exists for a broad range of 
pump values, and persists until the pitchfork bifurcation, 
at P = 1. At this point, two new branches of Too bifur- 
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cations emerge. The upper branch (b) correspond to a 
homocUnic bifurcation connecting one saddle with itself, 
while the lower branch (c) corresponds to a heretoclinic 
connection between two symmetric saddles. Three phase 
portraits corresponding to the different global bifurca- 
tions are shown in Fig. 2. The simultaneous presence 
of the three reported types of global bifurcations is a un- 
usual phenomenon. We note that a similar picture results 
from the analysis of two coupled van der Pol oscillators 
with delay coupling , however in a quite different con- 
text. 
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FIG. 2: Phase portraits of the acoustic field for parameters 
close to the different global bifurcations; (a) gluing, for i] — 
0.65 and V = 1.22, (b) homoclinic, for 77 = 0.7 and P = 0.96, 
and (c) heteroclinic, for rj = 0.65 and V — 1. 

The period of limit cycle solutions, as the pump param- 
eter is varied, is shown in Fig. 3 for two representative 
values of the nonlinearity parameter 77, below (dashed 
line) and above (continuous line) the codimension-three 
point where all the global bifurcations coalesce, at 77 ~ 
0.67. As expected from Fig. 1, at ?/ = 0.7 there is only 
one Too bifurcation, consisting in a homoclinic loop as 
shown in Fig. 2(b). For rj = 0.65, both gluing bifur- 
cation [Fig. 2(a)] and heterochnic connection [Fig. 2(c)] 
coexist. It is remarkable that this situation has been also 
described for a per iodically forced Navier-Stokes flow in 
hydrodynamics |1!tI ]. 

One prominent property found in some systems pre- 
senting homoclinic dynamics is excitability. As stated 
above, excitable systems present a highly nonlinear re- 
sponse to an external stimulus, with a well defined ex- 
citability threshold and a constancy in the reaction when 
perturbed above threshold. In addition, a refractory time 
has to elapse before the system can be excited again. 
These neuron-like properties can be also present in the 
magnetoacoustic resonator considered here, as we demon- 
strate below. 

The key for achieving excitability is the coexistence of 
a global bifurcation and a stable fixed point j3]. This 
occurs in the dark-shadowed region in Fig. 1. The ex- 
citability properties of a system can be characterized in 
several ways, in terms of its respose to different kinds 
of external perturbations. In order to demonstrate the 
main signatures described above, we consider first the 
behaviour of the system under a short (delta-like) per- 
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FIG. 3: Period of time-dependent solutions depending on the 
pump, for 77 = 0.65 (full line) and r] = 0.7 (dashed line). 
The period diverges to infinity at some values, denoting the 
presence of gluing (a), homoclinic (b) and heteroclinic (c) 
bifurcations. 

turbation. Figure 4(a) shows the amplitude response of 
the acoustic field after four perturbations with increas- 
ing amplitude, for parameter values i] = 0.7, V = 0.94 
and 7 = 0.1 (close to the homoclinic bifurcation bound- 
ary, curve (b) in Fig.l). The system, initially at rest, is 
perturbed at t = 0. The amplitude of the weakest pertur- 
bation was chosen to be below the excitability threshold 
yth (determined by the distance from the node to the sad- 
dle point), and the system relaxes smoothly to the rest. 
Three perturbations above the threshold, with different 
amplitudes, however generate three identical pulses or 
spikes. The amplitude of the perturbation only affects to 
the response time of the system: when stronger the stim- 
ulus, the smaller the time needed to develope a pulse. 

The delay time has been measured as the interval 
between the stimulus and the instant where the pulse 
reaches the maximum amplitude, for different overthresh- 
old perturbations. The results are shown in Fig. 4(b). 
The inset corresponds to the logarithmic representation 
of the numerical data, and demonstrates that the scal- 
ing law for the response time is ruled by the same law as 
the period of limit cycles close to homoclinic/heteroclinic 
bifurcations, i.e. 

T delay ^ -■J^^iV ~ yth) + C, (5) 

where A again corresponds to the unstable eigenvalue of 
the saddle point. For the parameters in Fig. 4, from 
the linear stability analysis results A = 0.01165, in good 
agreement with the value A — 0.01176 obtained from the 
slope of the linear fit in Fig. 4(b). 

Together with the excitation of single spikes, other 
characteristic features of many excitable systems are 
bursting and synchronization phenomena [2]. Burst- 
ing is the typical firing pattern displayed by neurons, 
and consists in the periodic emission of short trains of 



4 




FIG. 5: Bursting pattern induced by a periodic modulation 
of the control parameter. 

threshold (main signature of homoclinic dynamics) have 
been observed. 

The work was financially supported by the Spanish 
Ministerio de Educacion y Ciencia, and European Union 
FEDER (Project FIS2005-07931-C03-03). Discussions 
with Y. Fetisov are gratefully acknowledged. 

200 400 600 

delay time 




FIG. 4: Amplitude of the acoustic field in response to four ex- 
ternal perturbations of increasing amplitude, for i] = 0.7 and 
V = 0.94 (a). Delay time (dots) as a function of the strength 
of the perturbation. The inset corresponds to a logaritmic 
representation of the numerical data, and the continuous line 
to Eq. (6). 



fast spike oscillations, intercalated by quiescent intervals. 
Some systems (e.g. the neuron) present autonomous 
bursting, while in some others it can be induced by a 
weak periodic modulation of the control parameter, as 
has been shown in the CO2 laser with feedback 
We solve Eqs. (0J with a pumping term in the form 
V{t) = 7^(1 + mcos{ujpt)). Different bursting regimes 
have been observed depending on the modulation depth 
m and frequency ujp. Figure 5 shows a periodic burst- 
ing regime for V = 1, LUp = 0.005 and m = 0.0505. In 
this case, one burst is excited every period of the exter- 
nal driving (1:1 locking). Other modulation parameters 
result in different n : m phase locking regimes (defined 
by the ratio of bursting to modulation periods) , or even 
in chaotic bursting patterns, in agreement with other ex- 
ternally modulated excitable systems. 

In conclusion, we have presented for the first time an 
acoustic system displaying excitability. The system is 
formed by a magnetostrictive material excited by an os- 
cillating magnetic field. Some predictions of the theoret- 
ical model reported in this letter are in agreement with 
recent experimental results in a nonlinear magnetoacous- 
tic hematite (a-Fe203) resonator. In particular, bista- 
bility ans selfpulsing including oscillating solu- 
tions whose period strongly depends on the pump close to 
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